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I. INTRODUCTION 



It is now well established [y-|9| that a dissipative force operates in the dynamics of a fissioning nucleus as it descends 
from the saddle point to the scission. Since the fission process is expected to follow the mean field dynamics when the 
excitation energy of the nucleus is not too high, the physical processes that give rise to dissipation in fission are usually 
identified as one body mechanisms. The systematics of experimental mean kinetic energies of the fission fragments 
were fairly well reproduced ||l|,Hj9| using an one body dissipative force in classical dynamical calculations. However 
it has also been observed that the mean kinetic energy of the fission fragments is not very sensitive to the details 
of the dissipative forces. Though two body processes are less favoured compared to one body processes in the long 
mean free path dominated mean field regime, both one and two body dissipations in classical dynamical calculations 
have been found |3,|9| to describe the systematics of experimental mean kinetic energies. Nix and Sierk, on the other 
hand, suggested g-]^] in their analyses of mean fragment kinetic energy data that the effective dissipation is about 
four times weaker than that predicted by the wall plus window formula of one body dissipation. 

It has been pointed out |^-^ in recent years that prescission neutron and 7-multiplicities are sensitive to the details 
of the friction force in nuclear fission. In particular, it was shown Q that one body dissipation is preferred to two 
body viscosity in order to describe the neutron multiplicity data. It was further observed that a shape dependent 
dissipative force is essential for simultaneous reproduction of experimental values of prescission 7-multiplicities and 
fission probabilities. Specifically, it has been observed that the dissipative force in a compact configuration of the 
fissioning compound nucleus should be much smaller than that predicted by the wall formula of one body dissipation. 

The wall formula was developed by Blocki et al. in a simple classical picture of one body dissipation in which 
energy is transferred from the nuclear surface motion to the nucleon motion as a result of frequent collisions of the 
nucleons with the nuclear surface. According to the wall formula, the rate of energy dissipation is given as 

Ewpit) = pv h^da, (1.1) 



where h is the normal component of the surface velocity at the surface element da while the nuclear mass density and 
the average nucleon speed inside the nucleus are denoted by p and v respectively. The wall formula was also obtained 
from a formal theory of one body dissipation based on classical linear response techniques . 

The wall formula was originally derived for idealized systems employing a number of simplifying assumptions such 
as approximating the nuclear surface by a rigid wall and considering only adiabatic collective motions. The validity 
of these assumptions were scrutinised [|ll|,Q in the framework of random phase approximation (RPA) damping. It 
was shown JlT| , |l2| that in the limiting situation where the above assumptions are realized, RPA damping coincides 
with the wall formula. These works further brought out the importance of including the 2p2h states for realistic 
applications. 

We are however of the opinion that since the wall formula captures all the essential ingredie nts o f dissipation under 



mean field dynamics, it is still possible to improve upon the dissipation rate given by Eq.(l.l) by examining its 
various assumptions. One such assumption of the wall formula concerns the randomization of the particle motion. It 
is assumed Q that successive collisions of a nucleon with the one body potential give rise to a velocity distribution 
which is completely random. In other words, a complete mixing in the classical phase space of the particle motion is 
required. We have recently discussed p3| the implications of these assumptions and obtained a modification of the 
wall formula by relaxing the full randomization assumption in order to make it applicable to systems in which the 
mixing in phase space is partial. Considering chaotic particle trajectories which arise due to irregularity in the shape 
of the one body potential and which are responsible for irreversible energy transfer, we modified the wall formula in 
Ref. @ as, 

-BcwwF = ppv / h'^da, (1.2) 



where is a measure of chaos (chaoticity) in the single particle motion. The chaoticity fi is 1 for cavities of highly 
irregular shap es i n which most of the par ticle s move in a chaotic fashion, and one recovers the original wall formula 



given by Eq.(Ll) for such systems. Eq.(L2) becomes useful howev er fo r systems in which the particle motion is 



partially chaotic giving rise to a chaoticity less than 1. In fact, Eq.(1.2) gave a quite successful description |l4] of 



energy transferred to a gas from the of its container oscillating with small amplitudes. In what follows, we shall use 



the term "chaos wei ght ed wall formula" (CWWF) for Eq.(1.2) in order to distinguish it from the original wall formula 
(WF) given by Eq.(|lll). 

In the present work, we shall consider an ideal gas in a cavity which will be subjected to a time dependent quadrupole 
deformation. We shall calculate the energy transferred to the particles from the wall motion. We shall then compare 
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this transferred energy with that predicted by the chaos weighted wall formula (Eq.(p^)). In this model, the cavity 
will be progressivly deformed until the neck radius vanishes. The motion will be thus tailored to resemble that of a 
fissioning nucleus as it descends from saddle to scission. Our motivation in this work is twofold. First, we wish to test 
the applicability of the chaos weighted wall formula for large amplitude motions such as fission. Second, we would 
like to compare the chaos weighted wall formula with the phenomenological reduced friction suggested by Frobrich et 
al. 0. We shall find that the chaos weighted wall formula can account for the energy damping fairly accurately in 
large amplitude collective motions in an ideal gas. We shall also observe a qualiatative similarity between the chaos 
weighted wall formula and the phenomenological reduced friction. 

The dynamics of independent particles in time dependent cavities has been extensively studied earlier by Blocki 
and his coworkers |l^-[2^]. Considering classical particles in vibrating cavities of various shapes, a strong correlation 
between chaos in classical phase space and the efficiency of energy transfer from collective to intrinsic motion was 
numerically observed jist . In particular, it was demonstrated ||l^,|l6| that while the energy transfer is much smaller 
than the wall formula limit in a cavity undergoing quadrupole vibration, it reaches the wall formula limit for higher 
multipole vibrations. Similar conclusions were also reached ]r^pl[| when the particle motion was treated quantum 
mechanically, though the quantal energy transfers were found to be suppressed compared to the classical ones. In 
the present work, we shall demonstrate that the chaos weighted wall formula can describe the energy transferred to 
classical particles when a cavity is strongly deformed starting from spherical shape. 

The paper is organized as follows. We shall describe the fissionlike cavity motion in the next section. A procedure 
to obtain the irreversible energy transfer to the gas from the cavity motion will be outlined in Sec. III. Section IV 
will contain a comparison between this irreversible energy transfer and the chaos weighted wall formula prediction. 
A summary of our main results will be presented in Sec.V where we shall also draw our conclusions. 



II. A MODEL FOR FISSIONLIKE CAVITY MOTION 



We shall consider a gas of noninteracting particles in an axially symmetric cavity which is undergoing a quadrupole 
deformation. The initial velocity distribution for the particles will be assumed to be uniform upto a certain maximum 
value V similar to that of a completely degenerate Fermi gas. The deformation of the surface will be proportional to 
the Legendre polynomials P2 and the surface at any instant t will be defined as 



>!(i)P2(cos6i)}, 



(2.1) 



where 



a{t) = OLQ — a-i coswt. 



(2.2) 



Here ao and 02 are parameters which define the amplitude of deformation and lo denotes its frequency while \(t) 
is a volume conserving normalization factor. According to Eq.( p.2| ), the rate of change of deformation will be given 
as d = sinwi. The cavity is thus initially at rest with a deformation (ap ~ 0^2) which subsequently changes faster 
with time till the deformation reaches olq sX tot — ^. If we choose ao = 2.0 at which the neck radius becomes zero, 
the above cavity motion can serve as a model for fission dynamics during its transition from saddle to scission. By 
choosing a suitable value for 02, the cavity shape at the saddle point can also be defined. 

We can now evaluate Eq.(1.2) for a quadrupole deformed cavity and obtain 




where 
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CWWF 



a2?7W/i(a(i))/(Q;(t)) sin^ cut, 



(2.3) 



(2.4) 



and 



F{z) = 



[(1 + aP2iz)f + a2p^'(z)(l - z2)]i/2 



(2.5) 
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In the above equations, the adiabacity index rj is defined |Q as the ratio of the maximum speed of the surface tips 
to the maximum particle speed v and will be given by 

r,= ——, (2.6) 

and Eq denotes the total initial energy of the unperturbed gas. It may further be noted that for a ^ 1, Eq.(^) 
reduces to 

3 

— 2"^2'7'^Ai(a(i)) sin^wt, (2-7) 

CWWF 

which was used earlier for small amplitude shape oscillations. In the present work however we shall use Eq.( ^.3| ) 
since we shall be dealing here with large amplitude deformations. The essential steps to obtain Eq.(^) can be found 
in the Appendix. Subsequently the energy dissipated from the wall motion after a certain interval t will be obtained 
as 

^^a2V f\{a{t'))l{a{t'))sm''ujt'd{ujt'). (2.8) 




Eo 



CWWF 



In order to calculate the dissipated energy using Eq.(2.8), the chaoticity ^{a{t)) is required for all deformations a{t) 



through which the cavity evolves with time. The chaoticity for each deformation will be obtained [ [I3[ by considering 
particle trajectories in a static cavity with the same deformation and distinguishing between the regular and chaotic 
trajectories. Following the procedure outlined in Ref. [p2| , p^ , the chaoticity will be determined by uniformly 
sampling of the trajectories which originate from the cavity surface. Fig. |^ shows the calculated values of fi which will 
be subsequently employed to calculate the dissipated energy according to the chaos weighted wall for mula . It may 



also be noted that the original wall formula dissipation will be obtained by simply putting /i = 1 in Eq. ( p.q ) 



III. REVERSIBLE AND IRREVERSIBLE ENERGY TRANSFERS 

In the dynamical calculation of energy transferred to the particles in a cavity from its surface motion, the classical 
equation of motion of the particles is numerically solved. To this end, the initial position and the velocity vectors 
(lying within a Fermi sphere of radius v in velocity space) of a particle are chosen at random and the trajectory of 
the particle is follwed in time allowing reflections whenever it encounters the oscillating walls of the cavity . By 
considering a gas consisting of a large (typically several thousands) number of such noninteracting particles in the 
cavity, the relative energy transferred to the gas Efj-ans (t) /Eq is obtained upto an interval given hy ut = ^ at which 
the deformation becomes maximum. Fig. ^ shows the transferred energy calculated for different values of rj. The 
deformation parameters chosen are ao = a2 = 2.0 which correspond to an initial spherical shape and a maximum 
deformation of a = 2.0 at which the neck radius is zero. 

A part of the transferred energy Etrans{t) is reversible and is elastic in nature. This reversible part arises due 
to the symmetries present in the single particle Hamiltonian of the system and has been discussed by Blocki et al. 
101 in detail. For volume conserving systems, the reversible energy depends only on the deformation of the system 
and is independent of the speed of the deformation. It is of interest to note here that this property is quite distinct 
from that of a real gas. An intrinsic randomization of particle velocities is always assumed for an idea gas 
in a real container giving rise to the ideal gas laws according to which the energy of a gas does not change under 
adiabatic volume conserving deformations. The reversibility of the transferred energy in an ideal gas in a model 
container depends on the degree of nonintegrability of the single particle Hamiltonian and is thus fully reversible for 
an integrable system. This feature is illustrated in Fig. ^ where the relative transferred energy is shown over a full 
cycle of deformation for a cavity undergoing harmonic spheroidal deformation. The surface of the spheroidal cavity 
considered is defined by 

where 

c(t) -i?o(l + a(0), (3.2) 
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and 



a{t) = Ro{l + a{t))-^, (3.3) 

such that it undergoes volume conserving oscillations and particles inside such a cavity constitute an integrable system. 
The deformation will be given as a{t) = ao — ao cosLut, where we have used ao = 1.0 in our calculation. We find in 
Fig. 1^ that the transferred energy is fully reversible and is independent of the speed of deformation. 

For mixed (partially chaotic and partially regular) systems such as cavities undergoing quadrupole deformations, a 
part of the transferred energy is reversible while the rest is irreversible or dissipative in nature. The chaos weighted 



wall formula (Eq.(F2)) describes this irreversible part of the energy transfer. We shall therefore proceed to obtain 
the dissipative part from the total transferred energy as follows. 

Decomposing the transferred energy Etrans into its reversible Erev , and the dissipative Ediss , components, we write 

EtransiVi ") = Erev[a) + EdtssiVi (3-4) 

The incremental energy dissipation due to an incremental adiabaticity index will therefore be given as, 

AEd^ss{Al^,a) = AEtrans{Ar],a), (3.5) 

where 

AEdiss{Ari,a) = Edissiv + ^V^f^) - Ediss{ri,a), (3.6) 

and 

AEtrans{Ari, o) = Etransiv + Arj, a) - EtraiisiVi ")• (3-7) 



We shall obtain AEtransiAi], a) using Eq.(3.7) by first evaluating Etransiv ~^ a-^d EtransiVi from classical 

dynamical calculations. According to Eq.(3.5), this quantity will ba same as AEdissiArj, a) which shall be compared 



with the chaos weighted wall formula prediction in the next section. 

We should remark here that one can ideally obtain the full dissipated energy EdissiVi'^) from Eq.(3^) by first 
calculating the reversible part as 

Erev{a) = lim Etrans iV^C^)- (3-8) 

This scheme is difficult to implement numerically though. As the motion becomes slower (small 77 ), computation 
time increases prohibatively and numerical instability sets in. Further, it may be recalled that an interesting observa- 
tion was made in Ref. [ pT[ where it was shown that dynamical correlation between successive collisions of a particle 
with the cavity wall can substantially enhance the energy transfer rate beyond the wall formula limit for very slow 
motions of the wall. Thus we shall not obtain the desired reversible energy transfer from calculations at very slow 
wall motions. We have therefore used finite values of ry (within adiabatic limit) in the above incremental method for 
our purpose. 



IV. RESULTS 



We shall first consider a cavity which starts deforming from a spherical shape and reaches a maximum deformation 
of a = 2.0. The relative energy transferred to the gas for this system at different values of 77 is given in Fig. ||. 
Following the procedure outlined in the preceeding section, the incremental energy dissipation for various values of 
A77 are then obtained. This en ergy dissipation is also calculated from both the chaos weighted wall formula and the 
original wall formula using Eq.(2^). Fig. ^ shows the comparison. It is observed that the chaos weighted wall formula 



prediction for energy dissipation is very close to that from dynamical calculation. On the other hand, the original 
wall formula overestimates the energy dissipation by more than 50%. 

We shall next consider a cavity for which we set = 2.0 and 02 = 1.5. This gives rise to an initial deformation 
of a = 0.5 at t = 0. The results of dynamical calculations for energy transfer are shown in Fig. |[ The incremental 
energy dissipations are subsequently obtained and the comparison with predictions of the chaos weighted wall formula 
and the original wall formula is shown in Fig. ^. The chaos weighted wall formula prediction remains reasonably 
close to the dynamically calculated values for lower values of At] (0.02 and 0.04). For higher values of Arj (0.06 and 
0.08), the agreement becomes somewhat poorer at larger deformations. Nevertheless, the chaos weighted wall formula 
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provides a much better description of the dynamical results compared to the original wall formula for all the cases 
considered here. 

We shall now calculate the reduced friction coefRcient which corresponds to the dissipation according to the chaos 
weighted wall formula. The reduced friction coefficient f3 is defined as /? = 7 /M where 

E = jfRl (4.1) 

Here q is half the distance between the centres of mass of the equal fragments on each side of the neck divided by 
Rq, M is the total mass of the gas and 7 is the friction coefficient. The value of q varies from q = 0.375 for spherical 
configuration to q = 1 .0 for scission. The friction coefficient 7 and hence the reduced friction coefRcient (3 can be 
evaluated from Eq.( |2.§| ) for both the chaos weighted wall formula and the original wall formula (/i = 1). The reduced 
friction coefficients thus obtained are plotted in Fig. 0. The phenomenological reduced friction coefRcient introduced 
by Frobrich ||,|^ is also shown in this figure. 

We note in Fig. ^that the friction strength is strongly suppressed in chaos weighted wall formula for spherical or near 
spherical configurations. This reduction in the dissipative force at the early stages of nuclear fission is an important 
consequence of the chaos weighted wall formula. It is of interest to note that a reduction of the wall plus window 
formula of one body dissipation was also suggested |^-^ earlier. Further, such a reduction was phenomenologically 
found to be essential in order to reproduce the excitation functions of both the prescission neutron multiplicities 
and the fission probabilities. The subsequent increase of the reduced friction coefficient (due to chaos weighted wall 
formula) at larger values of q essentially reflects the increase of chaoticity with deformation. Since chaoticity depends 
sensitively on the nonintegrability of the system, the reduced friction coefficient at large q in turn depends on the 
details of the shape considered. Thus the comparison of our present calculation for a model shape evolution with the 
phenomenological friction is a qualitative one, both showing an increase of the reduced friction at higher deformations 
starting with very small values at the spherical configuration. We further note that the phenomenological friction near 
scission is much stronger than the original wall or its present chaos weighted version (both merge near scission) . This 
can possibly be understood |^ if other effects such as the exchange mechanism ("window" formula) and temperature 
dependence of nuclear friction are considered. 

V. SUMMARY AND CONCLUSIONS 

In the preceding sections, we have applied the chaos weighted wall formula to an ideal gas in a cavity undergoing 
fissionlike shape evolution in order to test the validity of the former for large amplitude motions. To this end, we have 
dynamically calculated the energy transferred from the wall to the gas and compared its irreversible part with the 
prediction of the chaos weighted wall formula. In order to simulate the macroscopic dynamics of symmetric fission in 
our model, we have chosen a cavity initially at rest with an initial deformation which is then progressivly deformed 
till the neck radius vanishes. In all the cases considered above, we find a fair agreement between the results from 
dynamical calculations and the predictions of the chaos weighted wall formula. The original wall formula is found to 
considerably overestimate the energy transfer. 

We have also compared the reduced friction coefficient extracted from the chaos weighted wall formula with that 
prescribed |^] phenomenologically in order to reproduce the excitation functions of both the precission neutron multi- 
plicities and the fission probabilities. Both the friction coefficients are found to be qualitatively similar. An important 
feature of the phenomenological friction is its strong suppression for compact shapes which is also present in the chaos 
weighted wall formula. 

We can therefore conclude that the chaos weighted wall formula provides a fairly accurate description of one 
body dissipation in dynamical systems similar to fissioning nuclei. Further, the qualitative similarity between the 
phenomenological friction in nuclear fission and the chaos weighted wall formula provides evidence for one body 
nature of the dissipative force acting in a fissioning nucleus. 
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APPENDIX A: 



We shall briefly outline here the essential steps to arriv e at the wall formula dissipation rate for cavities with large 
quadrupole deformations which has been used in Eq.( |2.3| ). 
The wall formula is given as, 



E = pv ri^da, 



(Al) 

where n = normal component of the velocity of the surface element da at a point r — r{9,(j)). The wall of an 
azimuthally symmetric quadrupole cavity is defined by, 



Taking the surface velocity component along the normal, 

fi = r cos 9n , 

where 0„ is the angle between f and n and is given by. 



(A2) 



(A3) 



and where we have used. 



cosf 



ds [,2 + (^)2]l/2' 



Hence 



Now from Eq.(|A2|) we get 



and 



ds^de[r^ + {-) ] . 



E = 2npv 



[^2 + (^)2]l/2 



sin ede 



r = ^[aP2 - j{l + aP2) 



where P2iz) = dP2{z)/dz and z = cos0. Using the above in Eq.(A4), we obtain 



where 



F{z) 



E^'-^rE,z)dz 



[P2{z)-\{^){l + aP2{z))Y{l + aP2{z)y 



[(1 + aP2{z)y + a^P!,\z){l - z2)]i/2 
The above equations give the rate of energy transfer for large quadrupole deformations. 



(A4) 



(A5) 



(A6) 



(A7a) 



(A7b) 
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FIG. 1. Variation of the chaoticity with deformation for quadrupole shapes. Solid circles are calculated values and the line 
is to guide the eye. 

FIG. 2. Growth of relative energy transferred to the gas from the wall motion. Curves a, b, c, d and e in the lower panel are 
for 77 = 0.02, 0.04, 0.06, 0.08 and 0.10 respectively while curves /, g, h and i in the upper panel are for 77 = 0.03, 0.05, 0.07 and 
0.09 respectively. The cavity motion is defined by ao = 02 = 2.0. 

FIG. 3. Relative energy growth for a cavity undergoing spheroidal deformation showing the complete reversibility of the 
energy transferred. The plot represents calculated values for 77 = 0.02,0.04,0.06,0.08 and 0.10 which are all indistinguishable 
from each other. 

FIG. 4. Growth of incremental dissipated energy relative to Eq for different values of Ar]. The solid lines axe from dynamical 
calculation. The long dashed and short dashed lines are the original wall formula (WF) and chaos weighted wall formula 
(CWWF) predictions respectively. 



FIG. 5. Same as Fig.2 but for the cavity motion defined by ao = 2.0 and 02 = 1.5. 



FIG. 6. Same as Fig.4 but for the cavity motion defined by ao = 2.0 and 02 = 1.5. 



FIG. 7. Reduced friction coefiicients from chaos weighted wall formula (solid line), original wall formula (long dashed line) 
and from Ref.[7] (short dashed line). 
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